Magnetized Neutron Star Mergers and Gravitational Wave Signals 
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We investigate the influence of magnetic fields upon the dynamics of, and resulting gravitational 
waves from, a binary neutron star merger in full general relativity coupled to ideal magnetohy- 
drodynamics (MHD). We consider two merger scenarios, one where the stars have aligned poloidal 
magnetic fields magnetic fields and one without. Both mergers result in a strongly differentially 
rotating object. In comparison to the non-magnetized scenario, the aligned magnetic fields delay 
the full merger of the stars. During and after merger we observe phenomena driven by the mag- 
netic field, including Kelvin-Helmholtz instabilities in shear layers, winding of the field lines, and 
transition from poloidal to toroidal magnetic fields. These effects not only mediate the production 
of electromagnetic radiation, but also can have a strong influence on the gravitational waves. Thus, 
there are promising prospects for studying such systems with both types of waves. 
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I. Introduction: To the long list of spectacular elec- 
tromagnetic signals which researchers have employed to 
study our universe, a new counterpart of gravitational 
origin will soon be added. Gravitational waves will both 
complement our ability to probe a variety of sources and 
in many cases provide the only option. These waves will 
likely be detected in the near future by a new genera- 
tion of laser interferometers [H, H, 0] • Among their most 
promising sources are merging binaries composed of black 
holes or neutron stars. These waves carry information 
about the dynamic, strong gravitational field near the 
binary. When one member of the binary is a neutron 
star, these waves can have a richer structure and carry 
additional information about the matter. For example, 
analysis of such mergers can potentially yield constraints 
on the nuclear equation of state and, as we show here, key 
insights on the stars' magnetic fields. The combination of 
gravitational wave data with other astronomical observa- 
tions in the electromagnetic spectrum will provide an en- 
hanced understanding of many energetic, strongly gravi- 
tating phenomena, such as (short) gamma ray bursts. 

Precise theoretical models are required to extract in- 
formation from gravitational waves. Recently, significant 
progress has been made in modeling the radiation ex- 
pected from binary black holes, binary neutron stars, and 
black hole- neutron star binaries. However, much less is 
known about how magnetic fields may alter the dynamics 
of such mergers, whose effects would be encoded in the 
gravitational waveforms. Magnetic fields are present in 
many astrophysical systems, and are important in several 
interesting phenomena, such as AGNs, pulsars, gamma 
ray bursts, etc. Observational evidence indicates that 
neutron stars have some of the strongest known mag- 
netic fields, with estimated field strengths ranging from 
10^^ G to as large as 10^^ G in magnetars [4]. Magnetic 
field effects in merging/collapsing scenarios have been in- 



vestigated in related systems. For instance, perturbation 
theory has been employed to examine magnetized dust 
collapse [H, Newtonian gravity to study magnetized bi- 
nary neutron star mergers , Post-Newtonian (PN) cal- 
culations to inspect the pre-merger dynamics [71] and full 
general relativity (reduced to 2D by symmetries) to ana- 
lyze axisymmetric differentially rotating star collapse Q . 
We study the merger in full, 3D general relativity with 
a framework able to study the fluid's behavior and ex- 
tract the gravitational waves in the wave-zone. Clearly 
this system is highly relevant in the mechanism underly- 
ing short gamma ray burst phenomena. However, since 
we do not include radiation transport our focus in the 
present work is on the gravitational wave signature and 
the fluid's behavior in a scenario probing the most ex- 
treme possible influence of the magnetic field. 

Magnetic fields can have a strong influence on the dy- 
namics of stars by introducing a number of instabili- 
ties (see e.g. (sl). These instabilities allow for energy 
exchange between the fluid and magnetic modes and 
can redistribute the fluid's angular momentum through 
different mechanisms (e.g., magnetic winding and brak- 
ing, magneto-rotational instabilities), depending on the 
strength of the magnetic field and the typical dynamical 
time scales [lol, • These effects can have a non-trivial 
effect on the (time dependent) multipolar structure of the 
source, and consequently on the resulting gravitational 
waves. A detailed understanding of magnetic effects in 
binary mergers and stellar collapse will provide valuable 
information not only for deciphering the observed gravi- 
tational waves, but they might also influence the design 
and tuning of future advanced detectors lH, 

In the present work we investigate possible magnetic 
effects in general relativity by simulating a merger of two 
compact stars with and without the presence of a mag- 
netic field. We study the merger in full, 3D general rela- 
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tivity with a framework able to study the fluid's behavior 
and extract the gravitational waves in the wave-zone. 
II. Overview of the numerical approach: We model 
the neutron star material using relativistic ideal mag- 
netohydrodynamics (MHD), coupled to the full Einstein 
equations of general relativity to represent accurately the 
strong gravitational effects during the merger. Our nu- 
merical techniques for solving these coupled equations 
have been tested previously 3, El, El, [ill- Some essen- 
tial elements of our implementation include the following. 
We use a computational infrastructure that provides dis- 
tributed Adaptive Mesh Reflnement (AMR) [13, 
with full sub-cycling in time, with a novel treatment of 
artiflcial boundaries 0]. The reflnement regions are de- 
termined by the truncation error of the solution during 
the evolution. We apply reflned boundary conditions at 
the outer boundary to ensure consistency of the solution 
over the computational domain H, [5^ and adopt high 
resolution shock-capturing methods together with diver- 
gence cleaning to control the no-monopole constraint . 
Our methods to solve the Einstein equations 1^, 13] and 
MHD equations Q are described elsewhere. 

Physical set-up: To begin exploring possible mag- 
netic effects in binary systems, we compare two different 
binary neutron star mergers, one with a magnetic fleld 
and one without. All other initial data are identical. 
For the magnetized stars, we choose a particular case 
in which each star initially has a poloidal magnetic fleld 
anti-aligned with the orbital angular momentum. Be- 
cause the stars' magnetic flelds are themselves aligned, 
this conflguration induces a magnetic repulsion as the 
stars approach. Strong fleld repulsion should impact the 
dynamics of the system and our results are compared 
with PN results presented in [7]. 

We construct the initial data with two identical neu- 
tron stars with zero spin angular momentum. Each star 
has a mass of M = 0.89 M©, a radius of 16.26 km, and 
a central density of 3.24 x 10^^ g/cm^. The stars are ini- 
tially separated by 60 km and are given initial tangential 
velocities slightly below the Keplerian value. We adopt a 
simple poly tropic equation of state choosing F = 2 (poly- 
tropic index n = 1) to model stiff nuclear matter [l5|. 

The poloidal magnetic fleld is calculated from the vec- 
tor potential = tu^ max(P — Pvac,0), where P is 
the pressure, w the cylindrical radius, and Pyac/c^ — 
10^ gm/cm^ ^2^. The fleld is initially conflned to the 
stellar interior, and its maximum magnitude is approx- 
imately 9.6 X 10^^ G. While this value is large for as- 
trophysical stars (excepting magnetars), it is small com- 
pared to the gravitational binding energy of the stars. As 
the individual stars do not rotate differentially, the initial 
fleld remains essentially constant until the stars are close 
enough for strong interactions. 

The initial data are evolved in a cubical computational 
domain given by [—1540 km, 1540 km] in each direction 
with an AMR conflguration having 7 levels of reflnement 



for the flnest resolved simulation (5 and 6 levels of re- 
flnement were used for comparison). The highest reso- 
lution region follows each star and has a grid spacing of 
A = 0.46 km. We extract the gravitational wave infor- 
mation and compute the strain {h-^^hx}- Observers are 
placed at three different radii from the center of mass, 
namely r = {440, 590, 740} km. These correspond to 
r ^ {5, 7.5, lOjAmax, where Amax is the largest wave- 
length in the observed gravitational wave. The extrac- 
tion takes place in the wave-zone, which is evident by the 
excellent agreement among the computed waveforms. 
IV. Results: The evolution of both binary pairs pro- 
gresses quite similarly. The stars' orbits shrink as en- 
ergy is lost through gravitational wave emission. As they 
merge, a bar is formed, which itself is unstable. The 
bar expels matter and collapses to form a differentially 
rotating, hypermassive star. While rotation initially sup- 
ports the star, the continued loss of energy through wave 
emission causes it eventually to collapse, forming a black 
hole. Comparing these two mergers, however, we see that 
the magnetic fleld strongly influences the details and time 
scales of the collapse, resulting in a solution that is clearly 
distinguishable from the unmagnetized case. 

Initially, the two solutions are nearly identical. When 
the stars begin to disrupt and shed matter, initiating the 
merger phase, the aligned magnetic flelds of the stars 
repel each other, delaying the flnal merger by 1-2 rota- 
tion periods with respect to the unmagnetized case. This 
delay is in agreement with post-Newtonian results pre- 
sented in [3]. As the merger proceeds, a number of rnag- 
netically induced phenomena are observed [1, [13, [llj . 
First, a shear interface forms giving rise to Kelvin- 
Helmholtz instabilities which are responsible for a rapid 
growth in the magnetic fleld strength, as shear energy 
is converted into magnetic energy (see Fig. [T]). Second, 
differential rotation of the resulting merged star winds 
the magnetic fleld up, producing a toroidal fleld and lin- 
early amplifying it as the star rotates. This trend sat- 
urates on the order of an Alfven time scale, which in 
our case is — 20Pc, where Pc is the rotational period 
of the star's core. In our simulation this period is about 
Pc = 0.735 ms (cf. [11). Third, at least two types of Tayler 

One is responsible for 
H. The sec- 
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instabilities are displayed 
inducing an e*^"^, m = 1 perturbation 
ond one, induces a magnetic buoyancy instability which 
sheds some portions of the star in low density regions 
(about 1% of the star's maximum density) which mani- 
fest in cavities that can be observed in the magnetic case 
but are absent otherwise. Finally, as the star collapses 
to a black hole, the magnetic fleld strength grows due to 
the compression of the fleld lines. The effects mentioned 
above strongly impact the fluid's distribution. Contrast- 
ing the magnetized case with the unmagnetized. Figure [2] 
shows density contours while Figure [3] displays the ratio 
of the equatorial semi-minor and semi-major axes. 
The combined effects of the magnetic fleld on the 



merger solution have a strong influence on the gravita- 
tional waveforms. Figure [4] shows the / = 2, m = 2 mode 
of the gravitational wave signal. At early times be- 
fore merger, the stars' magnetic flelds do not strongly 
interact, and the waveforms are nearly identical. The 
delayed merger is clearly visible in the large amplitude of 
h-^ which extends to t :^ 12 ms, or 1-2 more cycles than 
the non-magnetized case. Magnetic effects in the post- 
merger waveform are also evident: /i+ has both a smaller 
amplitude and lower frequency, especially for t > 20 ms. 
The magnetic fleld transports angular momentum away 
from the center, and as a result the rotation profile of 
the magnetized case displays a lower/higher frequency 
at the inner /outer regions of the star and a rather uni- 
form rotation profile in between together with a more 
axisymmetric configuration (see Figs. [2] and [3|). Thus the 
core is more axisymmetric and has a longer rotation pe- 
riod in the magnetized case than in the non-magnetized 
scenario. It is important to note that while our initial 
configuration has a large magnetic field, magnetic fields 
have been observed to saturate near this value in merger 
simulations that begin with weaker fields {c^ lO^^G) [6|. 
The similarity in these post-merger field strengths and 
configurations leads us to believe that the dynamics will 
demonstrate analogous behavior. Finally, the frequency 
range of the obtained waveform is about 2.7-6.4 kHz. 
LIGO is sensitive to strains of 4 x 10~^^/a/Hz in this 
range [l[ . The waves shown in Figure [4] have amplitudes 
comparable to this noise at distances of 2 (4) Mpc for the 
non-magnetic (magnetic) cases. Advanced LIGO would 
be able to detect these sources ten times farther, includ- 
ing the Virgo cluster. 




FIG. 1: Fluid density isopycnic (at p = 6 X lO^^g/cm^) and 
magnetic field distribution (in a plane slightly above the equa- 
tor) at t = 4.4ms. Kevin-Helmholtz instability in the shear 
layer has produced circulation cells that amplify the toroidal 
magnetic field. The cavities at both trailing edges are due to 
magnetic pressure inducing buoyancy. 



V. Conclusion: We have presented a simulation of 
magnetized neutron stars with full 3D GR and MHD with 
sufficient resolution to capture key magnetically-driven 



3 




FIG. 2: Equatorial density contours and velocity field at two 
times for both the non- magnetized (left column) and magne- 
tized cases (right column). The top frames are at t = 4.9 ms, 
clearly shows the delayed merger caused by magnetic field re- 
pulsion. The bottom frames are at t — 18.6 ms. With no 
magnetic field, the final differentially rotating star displays a 
more bar- like structure (m=2), while the magnetized star is 
considerably more spherical with a slight off-set. 



phenomena. We observe the transfer of shear energy into 
magnetic energy and magnetic winding, which are re- 
sponsible for the creation and amplification of a toroidal 
field in the resulting hypermassive star. The strong dif- 
ferential rotation of this star transforms the poloidal field 
into a toroidal field, the former decreasing by an order of 
magnitude. The toroidal field redistributes angular mo- 
mentum, resulting in an inner core rotating more slowly 
than the non-magnetized case. While some of these ef- 
fects have been seen to varying degrees in recent models 
(2D-axisymmetric simulations in GR [llj, 3D SPH simu- 
lations without GR \d\ ) , our implementation allows us to 
observe the magnetic interactions during the merger, con- 
sistently accounting for the tightening of the orbit due to 
the emission of gravitational waves, and does not restrict 
magnetic effects due to symmetries. The merger delay re- 
ported here has a strong inffuence in the resulting wave- 
forms, since they have more relatively high-amplitude cy- 
cles. Additionally, the subsequent waveform displays ob- 
servable differences in amplitude and frequency due to 
different magnetic instabilities. In our results, the in- 
spiral (chirp) phase displays no significant difference but 
the merger and after- merger epochs do. In fact a cross- 
correlation between the two waveforms is equal to 0.98, 
0.63, and 0.59, when integrating from t = to t = 7.5, 
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FIG. 3: Equatorial axis-ratio (semi-minor/semi-major) vs. 
time for the non-magnetized (blue-dashed line) and mag- 
netized (black-solid line) measured at a density of 1.24 x 
10^^ g/cm^ . the results corresponding to the non-magnetized 
case has been shifted in time by t = 3.5 ms so both merg- 
ers coincide. Notice how the magnetized case approaches 
much quicker the unit value while the non-magnetized re- 
mains about c^i 0.78. This translates into stronger amplitude 
in the produced gravitational waves. 
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FIG. 4: Gravitational wave signal /i+ (/ = 2, m = 2 mode) 
for the magnetized (solid line) and unmagnetized (dashed) 
mergers. Before t ^ 7.5 ms, both signals are essentially the 
same. After t :^ 7 ms, however, the delayed merger of the 
magnetized stars results in a much stronger signal lasting un- 
til t ^ 12 ms. After merger, the redistribution of angular 
momentum by the magnetic field results in a more axisym- 
metric remnant whose core rotates more slowly. As a result, 
at late times the non-magnetized merger produces radiation 
with larger/higher amplitude/frequency than its magnetized 
counterpart. 



t = 12.5, and t = 25 ms, respectively. Searches for these 
types of magnetic effects in gravitational waves can give 



(upper-limit) population estimates of these and related 
binaries. Much more work will be required to map out 
the possible configuration space of field strengths and 
orientations, but the tantalizing differences observed in- 
dicate interesting results can be obtained and linked to 
gravitational and electromagnetic observations. 
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